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ABSTRACT 

Dnmt2 enzymes are conserved in eukaryotes, where 
they methylate C38 of tRNA-Asp with high activity. 
Here, the activity of one of the very few prokary- 
otic Dnmt2 homologs from Geobacter species (Gs- 
Dnmt2) was investigated. GsDnmt2 was observed to 
methylate tRNA-Asp from flies and mice. Unexpect- 
edly, it had only a weak activity toward its match- 
ing Geobacter tRNA-Asp, but methylated Geobac- 
ter tRNA-Glu with good activity. In agreement with 
this result, we show that tRNA-Glu is methylated in 
Geobacter M\ih'\\e the methylation is absent in tRNA- 
Asp. The activities of Dnmt2 enzymes from Homo 
sapiens, Drosophila melanogaster, Schizosaccha- 
romyces pombe and Dictyostelium discoideum for 
methylation of the Geobacter tRNA-Asp and tRNA- 
Glu were determined showing that all these Dnmt2s 
preferentially methylate tRNA-Asp. Hence, the Gs- 
Dnmt2 enzyme has a swapped transfer ribonucleic 
acid (tRNA) specificity. By comparing the different 
tRNAs, a characteristic sequence pattern was iden- 
tified in the variable loop of all preferred tRNA sub- 
strates. An exchange of two nucleotides in the vari- 
able loop of murine tRNA-Asp converted it to the cor- 
responding variable loop of tRNA-Glu and led to a 
strong reduction of GsDnmt2 activity. Interestingly, 
the same loss of activity was observed with human 
DNMT2, indicating that the variable loop functions 



as a specificity determinant in tRNA recognition of 
Dnmt2 enzymes. 

INTRODUCTION 

Several chemical modifications of the bases in deoxyri- 
bonucleic acid (DNA) and ribonucleic acid (RNA) have 
been discovered in the last decades of research, which 
include methylation, pseudouridylation, hydroxymethyla- 
tion, formylation and carboxylation (1-4). In eukaryotes, 
three families of enzymes with strong sequence similarity 
to DNA-(cytosine C5)-methyltransferases have been found, 
namely Dnmtl, Dnmt2 and Dnmt3 (5-7). Dnmt2 enzymes 
are highly conserved with homologs present from lower eu- 
karyotes like Schizosaccharomyces pombe to higher eukary- 
otes like Homo sapiens (8,9). These enzymes contain all mo- 
tifs characteristic for DNA methyltransferases, but lack the 
large N-terminal domain found in the Dnmtl and Dnmt3 
enzymes. Later, in a seminal paper by Goll et al , Dnmt2 en- 
zymes were reported to have a robust methylation activity 
on tRNA^^P at position C38 (6). After the initial discovery 
of tRNA^^P methylation by Dnmt2, many studies have con- 
firmed this activity in different organisms (10-16). However, 
even though Dnmt2 methylates transfer RNA (tRNA), it 
utilizes a DNA methyltransferase-like catalytic mechanism 
in the methyl transfer reactions (10), which also explains the 
conservation of the corresponding amino acid motifs in the 
catalytic center. 

Dnmt2 knockout studies have indicated a connection to 
stress response and methylation of small RNAs (13,17,18), 
non-random sister chromatid segregation in stem cells 
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(19), mobile element and RNA virus control in Drosophila 
melanogaster (20,21), paramutation in mice (22) and adap- 
tation to different growth conditions in lower eukaryotes 
(12,15,16). However, it is still unclear whether these effects 
stem from the loss of methylation of C38 in tRNA^^P or 
whether other Dnmt2 targets are responsible for the phe- 
notypes. Many of the studies have also investigated methy- 
lation of alternate tRNA substrates by Dnmt2. In Dic- 
tyostelium discoideum, Dnmt2 was reported to methylate 
tRNAGi^(CUC/UUC) and tRNAGiy(GCC) very weakly 
in vitro (15). The S. pombe Dnmt2 homolog, Pmtl, was 
shown to have a weaker methylation activity on tRNA^^^ 
as weh (16). Z). melanogaster Dnmt2 was found to methylate 
tRNA^^i(AAC) and tRNA^iy(GCC) in vivo in addition to 
tRNA^^P (13). However, in all enzymatic studies, tRNA^^P 
was the most preferred substrate for Dnmt2. The molecu- 
lar basis for the specific interaction of Dnmt2 with its tRNA 
substrates is not yet known. Previously, we have mapped the 
binding of tRNA'^^P to human DNMT2 by mutating several 
conserved lysine and arginine residues in human DNMT2 
and concluded that the anticodon stem/loop of tRNA'^^P is 
a main region of contact for human DNMT2 for the methyl 
transfer reaction by human DNMT2 (23). Supporting this 
notion, MuUer et al proposed a role of the C32, A37 and 
C40 nucleotides in tRNA recognition (15). 

A recent phylogenetic analysis revealed that the Dnmt2 
gene family was most likely derived from a prokaryotic 
DNA-(cytosine C5)-methyltransferase (9). As mentioned 
above, the enzyme is highly conserved in all eukaryotic 
phyla. Strikingly, there are only a few putative Dnmt2 ho- 
mologs in bacteria, one of which is found in Geobacter 
species (6,9,24), another putative representative is present 
in Holophaga foetida (Supplementary Figure SI). These en- 
zymes are clearly defined by the presence of a CFT amino 
acid motif (CFI in Holophaga), which is characteristic for 
Dnmt2 enzymes and involved in the tRNA interaction 
(9,10). A bioinformatics analysis indicated that the Geobac- 
ter Dnmt2 enzyme is likely the outcome of a horizontal gene 
transfer from eukaryotes to Geobacter (9), which is known 
for its abihty of lateral gene transfer (25,26). This conclu- 
sion is supported by the additional identification of a puta- 
tive Dnmt2 homolog in Holophaga, because these bacteria 
are not closely related to Geobacter {Holophaga is a member 
of the Fibrobacteres/Acidobacteria group while Geobacter 
belongs to the 8-Proteobacteria). 

As part of a larger initiative aiming to investigate the 
function of Dnmt2 in many species including humans, mice, 
flies and unicellular eukaryotes and to define its biologi- 
cal role, we set out to investigate the biochemical proper- 
ties and substrate recognition of Geobacter sulfurreducens 
Dnmt2 (GsDnmt2). Since it has been shown that eukary- 
otic Dnmt2 enzymes are specific for tRNA^^P C38 methy- 
lation, our work was started in this direction. However, 
we surprisingly found that GsDnmt2 specifically methylates 
Geobacter tRNA^^^ in vitro and in vivo and not Geobacter 
tRNA^^P. This represents the first example of a Dnmt2 en- 
zyme with a completely switched substrate specificity. Based 
on this finding, we further identify the variable loop as a 
sequence determinant in tRNA recognition of Geobacter 
Dnmt2 and the human enzyme. 



MATERIALS AND METHODS 

Cloning and site-directed mutagenesis 

The gene for G sulfurreducens, GsDnmt2 (GSU0227), was 
amphfied from genomic DNA of G. sulfurreducens obtained 
from DSMZ (Braunschweig, Germany) using a primer set 
flanking the gene (forward 5^-GCC GCA TAT GAG GGC 
GGT CGA GCT CTT CTG-3^ and reverse 5^-GAT CGG 
AAT TCT CAC CCC TCC TCA GCC GGT AAC-30. 
The amplified gene was cloned into pET28a(+) using Ndel 
and EcoRI sites and successful cloning vahdated by DNA 
sequencing. The catalytic site variant C74A of GsDnmt2 
was created by the megaprimer site-directed mutagenesis 
method as described previously (27). This residue corre- 
sponds to C79 in human DNMT2, which is a key catalytic 
residue, mutation of which has been shown to inactivate the 
enzyme (10). The primer used for the site-directed mutage- 
nesis was 5 -GTA AGG CTG CGC AGG AGG CGA C- 
3^ The presence of the desired mutation was confirmed by 
sequencing. Constructs for other Dnmt2s are as described: 
human DNMT2 (28), S. pombe Pmtl (16) and D. dis- 
coideum DnmA (15). The gene for D. melanogaster Dnmt2 
was amplified from complementary DNA (cDNA) using 
the following primer pair: 5^-GCG GTG GTG CTC GAG 
TTA TTT TAT CGT CAG-3^ and 5^-GCG GCA GCC ATA 
TGG TAT TTC GGG TCT TAG AA-3^ and cloned into 
pET28a+ using Xhol and Ndel sites. The cDNA was pre- 
pared using RNA isolated from 0-6-h-old embryos of D. 
melanogaster Berlin Wild strain. The cloned sequence was 
found to contain an insertion present in the Isoform A of 
the D. melanogaster Dnmt2 gene, which was subsequently 
excised by polymerase chain reaction (PCR) mutagenesis 
performed with the following primers: 5^-ATG CCC AAT 
TGG ATG GAC AAA TAG TTG CCG CCT TGG-3^ and 
5 -ACT AAA TAG TTC TAA GAC CCG AAA TAC CAT- 

y. 

Protein expression and purification 

The wild-type and mutant C74A GsDnmt2 were expressed 
as His6-tagged proteins in Escherichia coli Rosetta 2 (DE3) 
cells. Protein expression and purification was performed 
as described previously (28). Briefly, transformed Rosetta 
2 (DE3) cells were grown to OD600 = 0.6. Protein ex- 
pression was induced with 1-mM IPTG (Isopropyl-p-D- 
1-thiogalactopyranosid) and cells were further incubated 
at 37°C for 3 h. The cells were sonicated in lysis buffer 
(30-mM potassium phosphate pH 7.5, 500-mM KCl, 0.1- 
mM DTT (dithiothreitol), 10-mM imidazole and 10% glyc- 
erol). Recombinant proteins were purified on Ni-NTA 
beads (Genoxon) and eluted with 200-mM imidazole in 
lysis buffer. The protein was then dialyzed against dial- 
ysis buffer I (30-mM potassium phosphate pH 7.5, 300- 
mM KCl, 10% glycerol, 1-mM ethylenediaminetetraacetic 
acid (EDTA) and 0.1-mM DTT) and dialysis buffer II (30- 
mM potassium phosphate pH 7.5, 100-mM KCl, 60% glyc- 
erol, 1-mM EDTA and 0.1-mM DTT). A similar proce- 
dure was followed for human DNMT2 purification. For D. 
melanogaster Dnmt2, the protein expression was conducted 
at 28° C. For the DnmA recombinant protein, expression 
was carried out at 22° C. The amounts of all recombinant 
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proteins were determined by their specific absorbance and 
the concentrations and purities were verified by Coomassie- 
stained sodium dodecyl sulphate-polyacrylamide gel elec- 
trophoresis (SDS-PAGE) gels. 

In vitro transcription of tRNA substrates 

The tRNA sequences were derived from the Genomic 
tRNA data base (http://gtrnadb.ucsc.edu/) (29). The tRNA 
substrates for methylation reactions were prepared by in 
vitro transcription basically as described (10). The DNA 
templates for the in vitro transcription were purchased from 
MWG. The template DNA and primer sequences are given 
in Supplemental Table S2. The templates were amplified by 
PGR with T7 and tRNA-specific primers and the yield and 
purity of the PGR was inspected on a native 8% acrylamide- 
TPE (Tris-Phosphate-EDTA) gels. In vitro transcriptions 
were performed in a total volume of 200 |jl1 containing 50 
|jl1 of 4x transcription buffer (200-mM Tris-HGl pH 8.1, 
5-mM spermidine, 25-mM DTT, 0.05% TritonX-100, 50- 
mM MgGl2 and 0.2-mg/ml bovine serum albumin), 20-mM 
NTPs (nucleoside triphosphates) mix, 1 U/|jl1 of RNasin 
plus (Promega), I U/|jl1 of T7 RNA polymerase (Fermen- 
tas) and 35 |jl1 of PGR-amplified template DNA. The reac- 
tion was incubated at 37°G for at least 3 h, followed by addi- 
tion of 15 U DNase I and incubation for another 30 min at 
37°G. Then tRNAs were separated on a 12% acylamide/7- 
M urea gel and the tRNA bands were excised from the gel 
(Figure IB). The tRNA was eluted with elution buffer (50- 
mM Tris-HGl pH 7.5, 300-mM sodium acetate and 0.5% 
SDS) overnight at room temperature followed by ethanol 
precipitation. 

In vitro methylation kinetics 

The in vitro methylation reactions of tRNAs were carried 
out as described before (30) with few modifications. Be- 
fore the methylation reactions, the substrate tRNAs were 
refolded by heating to 65°G and slowly cooled to ambi- 
ent temperature in the presence of 2-mM MgGl2 and 2 
U/|xl of RNasin plus (Promega). The subsequent methy- 
lation reaction was performed with 0.5 |jlM of the tRNA 
substrates and 1 |jlM of the Dnmt2 enzyme in methylation 
buffer (20-mM Tris-HGl (pH 8), 20-mM NH4OAC, 2-mM 
MgGl2, 2-mM DTT and 0.02-mM EDTA). The reactions 
were started by addition of 0.76-|jlM [methyl ^H]-AdoMet 
(Perkin Elmer) and carried out at 22° G with DnmA, 30° G 
with Drosophila Dnmt2 and 37°G with all other enzymes. 
The reactions were stopped at different time points by re- 
moving 5-|jl1 aliquots and pipetting them into 500 |jl1 of 5% 
trichloroacetic acid (TGA). The samples were then spot- 
ted onto a DE81 anion exchange filter paper disc (0.20 mm 
thick and 1.5 cm in diameter) and the unbound radioac- 
tive AdoMet was washed away with 5% TGA and absolute 
ethanol. The discs were then air dried and the bound ra- 
dioactivity was measured with a Hidex 300 SL liquid scin- 
tillation counter (Hidex, Finland) using Rotiszint eco plus 
(Roth, Germany). The reaction progress curves of the sin- 
gle turnover kinetics were fitted using Equation (1), which 
is derived from a single exponential model, 

CMPtheo(0 = /(I - e-^^') (1) 



with k\ being the single turnover rate constant of RNA 
methylation [min~^] and /being the saturation level of com- 
plete methylation (CPM). To obtain the initial rate of RNA 
methylation in GMP/min (vq). Equation (1) was differenti- 
ated at f = 0, which results in Equation (2): 

vo = f y<ki (2) 

In many reactions with less active enzyme/substrates com- 
binations a clear linear initial phase was observed. Then, 
the initial rate constant of RNA methylation was derived 
directly by linear regression of the reaction progress curve. 

tRNA binding analysis 

The tRNA-binding analysis was performed using 3^-end- 
pGp-labeled tRNAs. For labeling, 1 ixg of the in vitro tran- 
scribed Gs-tRNA^'P and Gs-tRNA^^^ was incubated with 
5 |jl1 [5'-P^2] Cytidine 3',5'-bisphosphate (pGp) (Hartmann 
Analytic) and T4 RNA hgase (NEB) in 50 |jl1 ligation buffer 
(50-mM Tris-HGl pH 7.5, 10-mM MgGl2, 1-mM DTT and 
1-mM adenosine triphosphate) and 10% DMSO (dimethyl 
sulfoxide) overnight at 4°G. The labeled tRNAs were puri- 
fied from unincorporated nucleotides by using Micro-Bio- 
Spin 6 columns (BioRad). The purified tRNAs were re- 
folded and incubated with increasing Dnmt2 concentra- 
tions in methylation buffer. After 30 min the reactions were 
spotted onto a nitrocellulose membrane and washed three 
times with the reaction buffer. The membranes were ex- 
posed to X-ray films and the intensity of the spots was mea- 
sured using ImageJ software. Data were fitted using the Ex- 
cel solver module to an equation describing a bimolecular 
binding equilibrium. 

Total RNA isolation from G, sulfurreducens 

G. sulfurreducens was grown anaerobically (N2:G02 80:20) 
at 30°G in NBAF medium with acetate (15 mM) as the elec- 
tron donor and fumarate (40 mM) as the electron acceptor 
as previously described (31). Total RNA isolation from G. 
sulfurreducens was carried out with the MasterPure RNA 
purification kit (Epicentre) following the manufacturer's in- 
structions. The isolated total RNA was treated with DNase 
I to remove contaminant genomic DNA. RNA concentra- 
tions were determined by ultraviolet spectroscopy. 

tRNA bisulfite sequencing 

The tRNA bisulfite sequencing was performed as de- 
scribed previously (32). For this study, tRNA^'P, tRNA^^^, 
tRNA^^^"^ and tRNA^^^ were selected because they con- 
tain cytosine at position 38 in their sequence. Briefly, 3 |xg 
of total RNA was treated with bisulfite following the in- 
structions of the manufacturer (EpiTect Bisulfite Kit, Qi- 
agen). The bisulfite-treated RNA was later used to generate 
a cDNA by reverse transcription using a stem loop primer 
(5^-GTG AAG TGG TGT GGT GGA GTG GGG AAT 
TGA GTT GAG TGG TAA GAA AATG-3^ 5^-GTG AAG 
TGG TGT GGT GGA GTG GGG AAT TGA GTT GAG 
TGG TAA TAA TGGG-3^ and 5^-GTG AAG TGG TGT 
GGT GGA GTG GGG AAT TGA GTT GAG TGG TAA 
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Figure 1. G. sulfurreducens Dnmt2 methylates C38 of tRNA^^P. (A) Methylation oi D. melanogaster tRNA^^P and murine-tRNA^^P with the GsDnmt2 
enzyme. Control reactions using the catalytically inactive GsDnmt2 mutant C74A and Dm-tRNA^^P C38U mutant show no methylation. (B) Methylation 
of tRNA^^P from G. sulfurreducens with GsDnmt2 showing weak methylation. Control reactions with the C38U mutant of the tRNA^^P or the catalytically 
inactive C74A variant of GsDnmt2 showed no methylation signal. (C) Summary of the tRNA methylation data reporting average methylation rates and 
SEM taken from two to three independent experiments. 
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Figure 2. GsDnmt2 prefers methylation of G. sulfurreducens tRNA^^^ over tRNA^^P. (A) Examples of in vitro methylation of Geobacter tRNA^^P and 
tRNA^^^ by GsDnmt2. (B) Average methylation rates and SEM based on three repeats of the experiments. Panel (B) also includes results of control 
reactions with tRNA^^P and tRNA^^^ C38U variants and with the catalytically inactive C74A GsDnmt2 protein. The error bars indicate the SEM from 
at least three different experimental repeats. 



TAA ACAC -30 for tRNA^^P, tRNA^^^ and tRNA^^^, re- 
spectively. The reverse transcriptase-PCR product was fur- 
ther amphfied by PGR using a universal stem loop primer 
5^-CAC GAC ACC AGT TGA-3^ and deaminated tRNA- 
specific primers with the following sequence: 5^ -TAG TTA 
AGT TGG TTAT-3^ 5^- GTT TAG TGG TTA GGA-3^ and 
5^-TTA GTT TAG TGG GAGA-3^ for tRNA^^P, tRNA^^^ 
and tRNA^^\ respectively. The PGR products were cloned 
by TOPO-TA cloning (Invitrogen, Germany). Colony PGR 
was performed to confirm the presence of desired inserts. To 
determine the RNA methylation, 60 clones were sequenced 
for each tRNA (33). 



RESULTS 

GsDnmtl is an active enzyme capable of methylating C38 in 
tRNA^^P 

After successfully cloning the GsDnmt2 gene into the 
pET28a+ vector, expressing and purifying the protein (Sup- 
plementary Figure S3), its activity was tested using in 
vitro transcribed tRNAs (Supplementary Figure S4). It was 
known from many previous reports that Dnmt2 from all the 
eukaryotic organisms methylate tRNA^^P (6,10,12-16). To 
test whether the bacterial Dnmt2 also methylates tRNA^^P, 
in vitro transcribed tRNA^^P of D. melanogaster was used 
as the initial substrate. Since the tRNA methylation rate of 
human DNMT2 is in the range of few turnovers per hour 
(10) and the Dnmt2 enzymes from other species studied 
here showed similar or lower activity, it was not possible 
to conduct steady- state methylation experiments. Therefore, 
we followed the initial time course of tRNA methylation re- 
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Figure 3. tRNA^^^ but not tRNA^^P is methylated at C38 in G. sulfurre- 
ducens. The figure shows resuhs of RNA bisulfite sequencing of tRNA^^P, 
tRNA^i^, tRNA^i^-2 and tRNA^^^ from total RNA isolated from G. sul- 
furreducens. The cytosine residues present in the corresponding tRNA are 
indicated in the top row. Each subsequent row represents an independent 
clone that was sequenced. Black boxes indicate methylated cytosines and 
white boxes indicate unmethylated cytosines. The arrows highhght the C38 
positions. 



actions under single turnover conditions using excess of en- 
zyme. As shown in Figure lA and C, wild-type GsDnmt2 
methylated tRNA^^^ while no activity was detected with the 
active site mutant (C74A) in which the catalytic cysteine 74 
residue, which is located in a highly conserved PCQ mo- 
tif in Dnmts, was altered to alanine. The corresponding ex- 
change in human DNMT2 has been shown to inactivate the 
enzyme as well (10). The observation that the C74A variant 
lost the methylation activity indicates that GsDnmt2 is an 
active tRNA methyltransferase and capable of methylating 
in vitro transcribed D. melanogaster tRNA^^P . The next aim 
was to confirm that the site of methylation by GsDnmt2 is 
the same as by other members of the Dnmt2 family. For this, 
a C38U mutant of the D. melanogaster tRNA^^P was pre- 
pared and methylation reactions were performed with the 
wild-type GsDnmt2 enzyme. A complete loss of methyla- 
tion was observed with this mutant tRNA, indicating that 
cytosine 38 is the methylation target site for GsDnmt2. Fi- 
nally, methylation of murine tRNA^^P by GsDnmt2 was 
demonstrated as well. In summary, these data indicate that 
GsDnmt2 methylates tRNA^^P from D. melanogaster and 
M mus cuius at C38. 



GsDnmt2 is weakly active on Geo bact er-encoded tRNA^^P 

After confirmation of the activity of GsDnmt2 on tRNA^^P 
derived from flies and mice, the catalytic activity of this 
enzyme on Geobacter-Qncodcd tRNA^^P was investigated, 
which shares large similarity with eukaryotic tRNA^^P 
(Supplementary Figure S2). For this, G. sulfurreducens 
tRNA^^P (Gs-tRNA^^P) was prepared by in vitro transcrip- 



tion and an in vitro methylation was performed using the 
recombinantly expressed wild-type GsDnmt2. However, as 
shown in Figure IB and C, GsDnmt2 was only weakly ac- 
tive on its cognate tRNA^^P. No activity was observed with 
the C74A mutant protein, confirming that the radioactive 
signal observed was due to tRNA methylation by wild-type 
GsDnmt2. In addition, the site of methylation was also con- 
firmed to be cytosine 38 in Gs-tRNA^'P by C38U muta- 
tion. The weak methylation of Gs-tRNA^^P was not re- 
lated to the preparation of the tRNA, because the same 
tRNA preparation was efficiently methylated by Dnmt2 en- 
zymes from other species (see below. Figure 4). In summary, 
when comparing methylation of tRNA'^^P from flies, mice 
and Geobacter, GsDnmt2 unexpectedly showed the weak- 
est methylation activity with its own Geobacter-QncodQd 
tRNA^^P. 



GsDnmt2 shows unusual preference for tRNA^*" over 
tRNA^^P 

After the initial observation that G sulfurreducens Dnmt2 
showed lower methylation of its Geobacter-Qncodod 
tRNA^^P than of tRNA^^P from Drosophila or M mus- 
culus, the question arose whether any other Geobacter 
tRNA might be more preferred as a substrate. To this end, 
Geobacter tRNAs were inspected for the presence of a cyto- 
sine at position 38 of the tRNA sequence. Using the UCSC 
Genomic tRNA Data Base (http://gtrnadb.ucsc.edu/) (29) 
we identified tRNA^^^ as well as two isodecoders of each 
tRNA^i% tRNA^i^ and tRNA^^^ as Dnmt2 candidate 
substrates in Geobacter, in addition to tRNA'^^P. Among 
them, both isodecoders of tRNA^^^ were selected for 
methylation analysis, because they show high similarity 
to tRNA^^P, as well as tRNA^^^ and one of the tRNA^^^ 
and tRNA^^^ sequences (Supplementary Table SI). All of 
these tRNAs were synthesized by in vitro transcription and 
assayed for methylation by GsDnmt2. Our results showed 
that tRNA^^P was weakly methylated as stated above 
(Figure 2), but the initial reaction rate for methylation of 
tRNA^i^ (tRNA32) was 10-fold higher than methylation 
of tRNA^^P. The methylation of the second isodecoder of 
tRNA^i^-2 (tRNA20) was intermediate between tRNA^'P 
and the other tRNA^^^ (Figure 2) and methylation of other 
tRNAs was not detectable (data not shown). These results 
motivated us to continue our study using the tRNA^^^ 
(tRNA32) isodecoder. A mutant tRNA^^^ carrying a C38U 
mutation was prepared and we observed that the mutation 
caused a complete loss of methylation, indicating that the 
methylation in tRNA^^^ happens only at position C38 
(Figure 2). In addition, the C74A mutation in GsDnmt2 
led to loss of methylation on tRNA^^^ (Figure 2). The 
preferential methylation of Gs-tRNA^^^ was also con- 
firmed using a gel-based methylation assay (10) in which 
the methylated tRNA is directly observed (Supplementary 
Figure S5). 

To test if the substrate preferences of human and Geobac- 
ter Dnmt2 were based on preferential binding of one type 
of tRNA, we determined their equilibrium constants of 
binding to Gs-tRNA^^^ and Gs-tRNA^'P (Supplementary 
Figure S6). Our data show that the human enzyme binds 
tRNA about four to five times stronger than the Geobac- 
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Figure 4. Dnmt2 enzymes from various species prefer metliylation of G. sulfurreducens tRNA^^P. (A) Examples of metliylation ]<:inetics are sliown for 
wild-type Dnmt2 enzymes from S. pombe (Pmtl), D. discoideum (DnmA), D. melanogaster and H. sapiens using in vitro transcribed Gs-tRNA^^P and 
Gs-tRNA^^^ and their respective C38U variants. (B) Average metliylation rates and SEM based on two repeats of the experiments. 



ter enzyme, which is in agreement with its general higher 
enzymatic activity (compare Figures 2B and 4B). How- 
ever, the binding of the two different substrates was simi- 
lar in both cases, indicating that ground-state tRNA bind- 
ing does not contribute to specificity of tRNA methylation 
by Dnmt2. A similar observation has been made previously 
in a study aiming to identify amino acid residues involved 
in tRNA binding and recognition of human DNMT2 (23). 
We conclude that GsDnmt2 shows a striking change in 
substrate specificity when compared with other Dnmt2 en- 
zymes, which is not based on tRNA-binding preferences. 

Cellular methylation of tRNAs in G. sulfurreducens 

As described above, we have found that GsDnmt2 has a dif- 
ferent specificity than other Dnmt2 enzymes in vitro. There- 
fore, the cellular methylation patterns of tRNAs in G. sul- 
furreducens were investigated. For this analysis, total RNA 
was isolated from wild-type G. sulfurreducens and the cy- 
tosine methylation investigated by bisulfite conversion, re- 
verse transcription and cloning of the converted tRNA^^P, 
tRNA^^^ and tRNA^^^ followed by sequencing of individ- 
ual clones. As shown in Figure 3, the methylation at the 
C38 position of tRNA^^^ in wild-type G. sulfurreducens was 
22%. The tRNA^^^-2 isodecoder was methylated by 36% at 
the same site. In contrast, signals for tRNA^^P and tRNA^^^ 
C38 methylation were at background levels of about 3%. 
Hence, both tRNA^^^ isodecoders were efficiently methy- 
lated in vitro and they are also methylated in cells, while 
tRNA^^P, which showed the weakest in vitro methylation 
among all substrates tested, was not methylated in cells. 



We conclude that cellular tRNA methylation at C38 mir- 
rors the in vitro specificity GsDnmt2. The same observa- 
tion has been made with human DNMT2. In human cells 
tRNA^^P but not tRNA^^^ is methylated at C38 (34,35), 
which is in agreement with the in vitro specificity of the hu- 
man DNMT2 enzyme. The incomplete methylation of the 
tRNA^^^ in Geobacter is similar to what has been observed 
for tRNA^'P in other species by RNA bisulfite (16,32,34). 
It is possible that the reverse transcription step of the RNA 
bisulfite method is more efficient with unmodified, imma- 
ture tRNAs, which may influence the results. 

Comparison of methylation specificities of Dnmt2 homologs 
from various species 

We concluded in the last paragraphs that GsDnmt2 strongly 
prefers to methylate Gs-tRNA^^^ not Gs-tRNA^'P both in 
vitro and in vivo. However, as shown above, tRNA'^^P from 
other species like D. melanogaster and M musculus were 
good GsDnmt2 substrates. This observation suggested that 
there may be changes in the Geobacter tRNAs that increase 
methylation of Gs-tRNA^^^ but reduce the reaction with 
Gs-tRNA^^P. Hence, the specificities of Dnmt2 enzymes 
from other species were investigated for the methylation 
of Geobacter tRNA^'P and tRNA^^^ also using the C38U 
mutant tRNA versions as controls (Figure 4). Dnmt2 en- 
zymes from different species that represent key eukaryotic 
phyla were expressed and purified (Supplementary Figure 
S3) and their methylation activities were analyzed. The as- 
says were done in parallel for a specific enzyme with all 
tRNA substrates. The results showed that human DNMT2, 
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Figure 5. Comparison of the specificities of Dnmt2 homologs from differ- 
ent species for methylation of Geobacter tRNA^^P and tRNA^^^. Data are 
replotted from Figures 2B and 4B and the relative activities of the Asp/Glu 
pairs were always normalized to the more active substrate for a better com- 
parison. Error bars were based on normalized individual experiments and 
correspond to the SEM of two to three repeats. 

S. pombe Pmtl and D. discoideum DnmA, all preferred the 
Gs-tRNA^^P substrate (Figure 4). D. melanogaster Dnmt2 
also showed some activity on Gs-tRNA^^^, but it still had a 
clear preference for Gs-tRNA^^P (Figure 4). The C38U mu- 
tations completely abrogated the methylation activity for all 
the enzymes, confirming that methylation occurred at posi- 
tion C38. Although the results with the other Dnmt2 en- 
zymes cannot be directly translated to their biological role, 
because the methylation was tested on G sulfurreducens tR- 
NAs, our data clearly demonstrate that the change in the 
biological target observed with GsDnmt2 is specific for this 
enzyme among all Dnmt2 enzymes tested here (Figure 5). 
These results show that the swap in specificity of GsD- 
nmt2 from Gs-tRNA^'P to Gs-tRNA^^^ must be due to 
changes within the GsDnmt2 protein and the Gs-tRNAs 
which affect the interaction of GsDnmt2 with Geobacter- 
encoded tRNAs. However, these alterations are only effi- 
cient when combined, because GsDnmt2 still methylates 
tRNA^^P from other species and Dnmt2 enzymes from 
other species prefer methylating tRNA^^P from Geobacter. 

The variable loop functions as a sequence determinant in the 
tRNA recognition of Dnmt2 

To investigate the molecular reason for the inversion of 
the substrate preference of the GsDnmt2 enzyme with its 
matching Geobacter tRNAs, we compared the secondary 
structures of tRNA^^P and tRNA^^^ from several species 
and noticed that Gs-tRNA^^P has a longer variable loop 
with an extra guanine (Figure 6). Furthermore, the variable 
loops of all preferred Dnmt2 substrates (namely tRNA^^P 
from all species but Geobacter and Gs-tRNA^^^) have at 
least one A at position 45 or 46 of the variable loop, while all 
the inactive counterparts contain a GG dinucleotide. Based 
on this, we speculated that the variable loop may be a speci- 
ficity determinant of GsDnmt2. To test this hypothesis, the 
variable loop of murine tRNA'^^P was swapped with that of 
murine tRNA^^^ by an exchange of two nucleotides (AG AC 



to CGGC) which introduced the GG dinucleotide sequence 
into tRNA^^P. The mutant tRNA^^P was synthesized by in 
vitro transcription. Equal amounts of the wild- type and mu- 
tant tRNAs were used in methylation reactions with GsD- 
nmt2 (Supplementary Figure S4). Strikingly, the alteration 
of these two bases in the variable loop of murine tRNA^^P 
led to a strong decrease in the methylation rate (Figure 7A). 
This result supports the hypothesis of an involvement of the 
variable loop in the enzyme-substrate recognition of GsD- 
nmt2. To further extend this finding, methylation reactions 
were performed with the same substrates using the human 
DNMT2 enzyme. Here, also, a similar reduction of activ- 
ity was observed with the mutant tRNA^^P compared to 
the wild type tRNA (Figure 7A), indicating that the vari- 
able loop is also important for tRNA^^P recognition by hu- 
man DNMT2. Next, we investigated if the variable loop of 
Gs-tRNA^^^ also has a role in recognition. As shown in 
Figure 7B, replacement of the original Gs-tRNA^^^ vari- 
able loop by either Gs-tRNA^^P or murine tRNA^^^ re- 
sulted in a strong reduction of its methylation by GsDnmt2. 
These data indicate that the GG dinucleotide in the variable 
loop functions as a universal anti-determinant for methy- 
lation by Dnmt2. However, replacing the variable loop of 
murine tRNA^^^ with the corresponding loop from murine 
tRNA^^P did not result in a significant increase of its methy- 
lation by human DNMT2 (data not shown), indicating that 
the variable loop is not the only recognition determinant of 
Dnmt2. 

DISCUSSION 

Dnmt2 enzymes are conserved in most eukaryotic species 
(6,9) and methylate tRNA^^P with high activity. Here, the 
activity of one of the few known bacterial putative Dnmt2 
enzymes was studied, which is found in Geobacter strains 
including G. sulfurreducens (GsDnmt2). We show that Gs- 
Dnmt2 actively methylates tRNA^^P from M. musculus and 
Drosophila that were used as model substrates in previous 
work. This result confirms that GsDnmt2 is a member of 
the Dnmt2 family as expected on the basis of the amino acid 
alignment. However, the endogenous Geobacter tRNA^^P 
was methylated only very inefficiently. Instead, GsDnmt2 
methylates Geobacter tRNA^^^ at the corresponding site, 
indicating a complete swap in its specificity (Figure 5) — 
a surprising and unexpected result that was confirmed in 
a methylation analysis of tRNA isolated from G. sulfurre- 
ducens. There are two isodecoders of tRNA^^^ in Geobacter 
(tRNA20, here called tRNA^^^-^, and tRNA32, here called 
tRNA^^^), which both were methylated in vitro and in vivo, 
although the relative methylation levels of the two isode- 
coders varied in the two data sets. In vitro, tRNA^^^ was 
methylated three times faster while in cells the methylation 
level of tRNA^^^"^ was 1 .6 times higher. There could be sev- 
eral explanations for this difference: (i) tRNA^^^"^ methy- 
lation could be less efficient under in vitro conditions, (ii) 
additional modifications could modulate the tRNA^^^ or 
tRNA^^^"^ methylation in vivo or (iii) additional modifica- 
tions could affect the recovery of methylated tRNA^^^ and 
tRNA^^^"^ differently after bisulfite conversion. However, 
in all data sets both isodecoders of Gs-tRNA^^^ were more 
methylated than Gs-tRNA^'P. 
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Figure 6. Clover leaf structure of the tRNA^^P and tRNA^^^ from different species. The target C38 is printed in bold and highlighted by a gray circle for 
the preferred substrate tRNA and by a white circle for the non-preferred one. Bases in the variable loop are colored. The sequences of human tRNA^^P 
and tRNA^i^ are identical to that of murine- tRNA^sp and tRNA^^^. 
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Figure 7. Methylation of tRNA variable loop mutants by GsDnmt2 and human DNMT2. (A) Methylation of murine-tRNA^^P with its original variable 
loop or with the variable loop of murine-tRNA^^^ by Geobacter Dnmt2 (left pair of bars) or human DNMT2 (right pair of bars). (B) GsDnmt2 methylation 
ofGs-tRNA^i^ with its original variable loop or with variable loops from Gs-tRNA^^P or murine- tRNA^i^.Thefi gure shows averages of the methylation 
rates normalized to the wild-type tRNA. Error bars indicate SEM derived from two repeats of the experiments. 



We have shown here that Geobacter Dnmt2 has an ~ 100- 
fold swapped specificity for Geobacter tRNAs, because it 
prefers tRNA^^^ roughly 10-fold over tRNA^'P, while all 
other Dnmt2 enzymes showed an at least 10-fold prefer- 
ence for tRNA^'P over tRNA^^^ (Figure 5). What could 
be the molecular basis of the unexpected loss of Geobac- 
ter tRNA^^P methylation by GsDnmt2 and the correspond- 
ing gain of methylation of tRNA^^^? So far, various amino 
acids involved in tRNA recognition of the human enzyme 
could be mapped (23), but the sequence elements needed 
on the tRNA side were largely unknown. Initially, the G34 
base at the wobble position of the anticodon, which is mod- 
ified to mannosylqueuosine in eukaryotic cells, was asso- 
ciated with Dnmt2 activity (6). However, the results with 
G. sulfurreducens are not in favor of that model, because 
Gs-tRNA^^P contains the G34, but it is a very weak sub- 
strate, while Gs-tRNA^^^, which is the preferred substrate. 



does not have a G34. By comparing substrate and non- 
substrate tRNAs of GsDnmt2, it became apparent that all 
non-substrates contain a GG dinucleotide in the variable 
loop as a characteristic feature. Indeed, introduction of this 
sequence into murine tRNA^^P and Geobacter tRNA^^^ 
drastically reduced their methylation by GsDnmt2. We can- 
not rule out the possibility that the two point exchanges in 
the variable loop might disrupt tRNA folding. However, 
the initial discovery that the GG dinucleotide is inhibit- 
ing GsDnmt2 was made with natural tRNAs, which are 
all fully functional in protein biosynthesis, which rules out 
loss of structure in these cases. Therefore, when taken to- 
gether, our results strongly suggest that the variable loop 
functions as an important specificity determinant in GsD- 
nmt2. Interestingly, the exchange in the murine tRNA^^P 
also reduced its methylation by the human DNMT2, indi- 
cating that the variable loop is a critical specificity determi- 
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nant for the tRNA recognition of human DNMT2 as well. 
However, introduction of the favorable variable loop into 
murine tRNA^^^ did not increase its methylation. Hence, 
our data indicate that Dnmt2 enzymes use further tRNA- 
specific sequence determinants in addition to the variable 
loop. These (so far unknown) contacts mediate the prefer- 
ence of the human DNMT2 enzyme for Gs-tRNA^^P, de- 
spite the presence of the unfavorable variable loop. Further- 
more, they prevent methylation of murine tRNA^^^, even 
after introducing a favorable variable loop. The contribu- 
tion of additional specificity determinants (besides the GG 
dinucleotide in the variable loop) to tRNA recognition by 
Dnmt2 enzymes is also supported by the observation that 
Dnmt2 enzymes from other species still strongly prefer Gs- 
tRNA^^P over Gs-tRNA^^^. 

The role of the variable loop as a specificity determi- 
nant in tRNA recognition by Dnmt2 is not without prece- 
dence. Examples for tRNA-interacting enzymes, which 
use the variable loop for tRNA recognition, include a 
bacterial tRNA (m7G46) methyltransferase (36) and a 
prokaryotic tRNA-dependent amidotransferase (37), and 
the same strategy is also used by several aminoacyl-tRNA- 
synthetases (38). Recognition of nucleic acids by proteins, 
in general, is based on direct and indirect readout, i.e. di- 
rect contacts of amino acids to nucleobases, which detect 
the nature of the base, and contacts to the nucleic acid back- 
bone, which detect the structure of the nucleic acid. The GG 
dinucleotide in the variable loop of the tRNA could inter- 
fere with either of these processes. However, the observa- 
tion that the variable loop of preferred tRNA substrates is 
only characterized by the absence of the GG dinucleotide 
but it does not contain a specific nucleotide sequence might 
argue against direct base contacts. This suggests that an in- 
direct readout mechanism is operational in which the pres- 
ence of the GG dinucleotide induces a structural change in 
the tRNA which then prevents its methylation. 

The finding that GsDnmt2 methylates tRNA^'P from 
other species (as all other members of the Dnmt2 family 
do) but not from Geobacter provides strong support for the 
hypothesis that the Gs-Dnmt2 was acquired by horizontal 
gene transfer, because there is no reason for this enzyme to 
methylate tRNA^^P, which is not a substrate in Geobacter 
itself In the horizontal gene transfer view, the activity to- 
ward tRNA'^^P can be interpreted as an evolutionary relict. 
Based on our data, one may speculate that initially a Dnmt2 
gene with tRNA^^P preference, but also the ability to methy- 
late tRNA^^^, was acquired by Geobacter. Precedence for 
such side activity can be seen in many species including 
D. melanogaster, S. pombe and D. discoideum (13,15,16). 
Methylation of the Gs-tRNA^^^ was favorable, which led 
to the stable integration of the Dnmt2, whereas methyla- 
tion of the endogenous tRNA^^P was lost in a co-evolution 
of GsDnmt2 and the Gs-tRNA^^P . At present, one can only 
speculate about the role of Dnmt2 in Geobacter. In previous 
studies, tRNA methylation by Dnmt2 was found to improve 
the stability of the tRNA (13) or affect tRNA charging (R. 
Shanmugam et al , in preparation), and prevent tRNA frag- 
mentation (13). Which of these roles, if any, is important for 
Geobacter tRNA^^^ cannot be said at present. The genera- 
tion of a Dnmt2 knock-out strain was unsuccessful (M. Ak- 
lujkar, unpubhshed information), which may indicate that 



Dnmt2 is an important gene in G. sulfurreducens, which is 
in line with the evolutionary arguments. This would be an 
interesting finding, since Dnmt2 enzymes are non-essential 
in other model organisms under normal growth conditions 
(6,12,13,15,21). 

CONCLUSIONS 

The first enzymatic characterization of the specificity of a 
bacterial Dnmt2 homolog, the enzyme from G sulfurre- 
ducens, revealed that it is the only known Dnmt2 enzyme 
that does not prefer methylation of tRNA'^^P. Instead, an 
unexpected swap of specificity of this enzyme was observed 
toward methylation of Geobacter tRNA^^^, which was con- 
firmed in vitro and in vivo. However, GsDnmt2 still methy- 
lates tRNA^^P from other species, and Dnmt2 enzymes 
from other species still methylate Geobacter tRNA^^P. We 
provide an evolutionary scenario that can explain this inter- 
esting finding. We noticed that the Geobacter tRNA^^P has 
a modified variable loop, which serves as a specificity deter- 
minant and causes the loss of methylation activity. Further- 
more, we show that this mechanism could be extended to 
the human enzyme as well, indicating that work with a bac- 
terial Dnmt2 homolog, finally, has led to important insight 
into the mechanism of the corresponding human enzyme. 
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